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Strong measurements give a better direct measurement of the quantum wave function 
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Weak measurements have thus far been considered instrumental in the so-called direct measure¬ 
ment of the quantum wavefunction [Nature (London) 474, 188 (2011)]. Here we show that direct 
measurement of the wavefunction can be obtained by using measurements of arbitrary strength. 
In particular, in the case of strong measurements, i.e. those in which the coupling between the 
system and the measuring apparatus is maximum, we compared the precision and the accuracy of 
the two methods, by showing that strong measurements outperform weak measurements in both for 
arbitrary quantum states in most cases. We also give the exact expression of the difference between 
the reconstructed and original wavefunctions obtained by the weak measurement approach: this will 
allow to define the range of applicability of such method. 


Introduction - In Quantum Mechanics the wavefunc¬ 
tion is the fundamental representation of any quantum 
system, and it offers the key tool for predicting the mea¬ 
surement outcomes of a physical apparatus. Its deter¬ 
mination is therefore of crucial importance in many ap¬ 
plications. In order to reconstruct the complete quan¬ 
tum wavefunction of a system, an indirect method, know 
as quantum state tomography (QST), has been devel¬ 
oped d]. QST is based on the measurement of comple¬ 
mentary variables of several copies of the same quantum 
system, followed on an estimation of the wavefunction 
that better reproduce the results obtained. This method, 
originally proposed for a two level system, has been ex¬ 
tended to a generic number of discrete quantum states [2] 
as well as to continuous variable state [3] . Recently Lun- 
deen et al. [3] proposed an alternative operational deh- 
nition of the wavefunction based on the weak measure¬ 
ment IBHH]. After the hrst demonstration, in which the 
transverse wavefunction of a photon has been measured, 
this method has been applied for the measurement of the 
photon polarization [9], its angular momentum jTO] and 
its trajectory mi. The method has been subsequently 
generalized to mixed states [12] to continuous variable 
systems [13] and compared to standard quantum state 
tomography in MM- 

By such method, that we call Direct-Weak- 
Tomography (DWT), a “direct measurement” of 
the quantum wavefunction is obtained: the term “direct 
measurement” refers to the property that a value 
proportional to the wavefunction appears straight from 
the measured probabilities without further complicated 
calculations or htting on the measurement outcomes [16] . 
As originally proposed [ 3 ], the method is based on the 
weak-measurement obtained by a “weak” interaction 
between the “pointer” (i.e. the measurement appa¬ 
ratus) and system. Weak measurements occur when 
the coupling between the pointer and the system is 
much less than the pointer width. As reported in the 
literature, “the crux of [the] method is that the first 
measurement is performed in a gentle way through weak 
measurement, so as not to invalidate the second or 
“Directly measuring [...] relies on the technique of weak 
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FIG. 1. Scheme of the original DWT method used to measure 
the wavefunction. 


measurement: extracting so little information from a 
single measurement that the state does not collapse”^. 

The interest about DWT is that the scheme in some 
cases may have experimental advantages over QST, in 
terms of simplicity, versatility, and directness [12] : it only 
requires a weak coupling of the system with an external 
pointer, a postselection of the hnal state of the system 
and a simple projective measurement of two complemen¬ 
tary observables of the pointer, a two-level system. QST, 
in contrast, requires measuring a complete set of noncom¬ 
muting observables of the system, which can be a very 
demanding requirement in systems with a large number 
of degrees of freedom. For instance the determination 
of the transverse spatial wavefunction of a single photon 
was hrst realized by DWT [ 3 ], as well as the measurement 
of a one-million-dimensional photonic state [5]. 

Here we show that the quantum wavefunction can be 
obtained by the same scheme used in DWT, but using 
only strong measurements: with this terms we here refer 
to measurements characterized by a strong coupling be¬ 
tween the system and the pointer. As explained below, 
a strong measurement does not always coincide with a 
projective measurement on the system. 

We thus demonstrate that the weak measurement is 
not necessary for the direct measurement of the wave- 
function. We then compare DWT with our method, 
showing that the use of strong measurements in most 
cases gives a better estimation of the quantum wavefunc¬ 
tion, outperforming DWT when both accuracy and preci- 
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sion are considered. Our analysis also allows to evaluate 
how the wavefunction estimated by DWT is related to 
the correct wavefunction, see eq. We also solved 

an unresolved question related to DWT: how “weak” the 
interaction should be such that DWT gives a correct esti¬ 
mation of the wavefunction. In particular, we will derive 
a sufRcient criterium for the applicability of DWT based 
on the measured probabilities, see eq. 

Review of Direet-Weak-Tomography - Let’s consider 
a d dimensional Hilbert space with basis {|a:)} with 
X — l,...,d. The states |x) are equivalent to position 
eigenstates of a discretized segment. A generic pure state 
in this basis can be written as 

d 

. ( 1 ) 

The scheme used in DWT is shown in hgure hrst, 
the following initial state |Tin) = l'0)x ® 
pared, with |0)p the pointer state. The pointer belongs 
to a bidimensional qubit space spanned by the states 
{|0)p) |l)p} [I3- The system is then evolved according 
to the following unitary operator: 

U^{e) = ^t-idTT^®ay, (2) 


where 9 is an arbitrary angle and tTx — |a:;)(cc|. The ap¬ 
proximation of the r.h.s. of eq. ([^ is obtained for small 0. 
The previous evolution corresponds to a pointer rotation 
conditioned to |'0)x being in the state |x). A projective 
measurement on the pointer, weakly coupled to the pho¬ 
ton position and followed by a projective measurement of 
the photon momentum allows to directly determine the 
wavefunction. Indeed, by post-selecting only the out¬ 
comes corresponding to the zero transverse momentum 
state Ipo) = Xlx I®)’ (unnormalized) pointer state 

becomes ^ ^[^|O)p+0^x|l)p] with ^ ^x- 

The choice of |po) is arbitrary, and a different value of 
the transverse momentum might be needed for particu¬ 
lar states, as explained below. Since a global phase is not 
observable, it is possible to arbitrarily clmose the phase 
of '0: we set the latter phase such that 0 is real valued 
and positive. In the hrst order in 9, the wavefunction can 
be derived directly as [1]: 


0m, a 


200 


(pW _ pw) p i(pm _ pw) 


)(a:) 


(x) 


)(®)' 


( 3 ) 


(x) 

where Pj represent the probabilities of measuring the 
pointer state into the diagonal basis |±)p = ;^(|0)±|1)), 
or the circular basis |L)p = l-^)p = 

■^(|0) — ^|1)). We note that, since the (real positive) 
proportionality constant is x independent, it can be 
obtained at the end of the procedure by normalizing the 
wavefunction. The different probabilities can be also ex¬ 
pressed in the framework of POVM, as detailed is SI. 


From now on, we indicate with 0m,x the (approximate) 
wavefurmtion obtained with the DWT method. We also 
dehne 0m = Y.x'^w^x = 
hx the global phase of 0m,x by 

Relation ^ was generalized to mixed states in [T^ . 
By repeating the measurements and changing the x pa¬ 
rameter in the evolution Ux{9), the full wavefunction can 
be reconstructed. We now show that a relation similar 
to 0 can be obtained by strong or arbitrary strength 
measurements. 

Arbitrary strength measurement - Measurement with 
arbitrary strength is obtained by choosing arbitrary value 
of 9 within 0 < 9 < tx/2. We start our analysis with 
strong measurements, corresponding to 0 = 7r/2. In 
this case the unitary operator ([^ becomes Ux{'xl2) — 
1 — |a:)(a;| (g) (1,^ + icTy). After the interaction, the initial 
state I Tin) is measured on the state \po)®\(f)f), where |0/) 
is the hnal polarization state. The amplitude for that 
transition is just A — (po|0)x(0/|O)-p — tfxVVdif/^fl-) V- 
This amplitude involves both the real and imaginary 
parts of 0a;, so its magnitude squared does too: by choos¬ 
ing different values of |0/), it is possible to determine the 
real and imaginary parts of 0x. In particular by choosing 
the hnal state |0/) as |l)-p, |+)-p, |—)-p, \L)'p and |R)-p 
states, the wavefunction can be obtained as: 


/ ^ 
0x = ^ 

20 


(^p{x) _ p{x 


2p(*) 


)+i(P 


(x) 

L 



(i) 

To obtain the above relation we hxed again 0 = |0|. 
It is very important to stress that, differently from the 
DWT method, the above result is exact, without any ap¬ 
proximation. We denote the previous relations as Direct- 
Strong-Tomography (DST) method. The difference with 
respect to the DWT is the need of measuring the pointer 
state also in the state ll)^. This extra requirement is 
compensated by the fact that the result is not approx¬ 
imated and the accuracy and precision of the method 
overcomes the DWT, as we will show in the following. 
We underline that the measurement in the \l)p state, 
and only in this state, corresponds to a projective mea¬ 
surement of the photon position, as the outcome of the 
measurement is proportional to |0xP (see S.I.). On the 
contrary, a projection of the pointer in the {|+),|—)} 
or {|L),|R)} bases acts as a partial quantum erasure 
on the which-position information: therefore a subse¬ 
quent momentum postselection allows to extract infor¬ 
mation about the real and imaginary part of 0^. As 
detailed in SI, for arbitrary 0, the wavefunction can be 
obtained as 3?e(0x) oc — P^^ + 2tan(|)p}*^ and 
Am(0x) oc - P^^\ 

Aeeuraey of DWT - In the case of DWT, the ob¬ 
tained wavefuction 0m,x is an approximation of the cor¬ 
rect wavefunction 0^. We now evaluate the accuracy of 
the DWT, namely the errors arising by using eq. ^ 
in place of the exact values of Q. As done in [T3], we 
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FIG. 2. Accuracy of the DWT: we show the probability pw of 
having ijjw < 0 and the probability px) of having an error T> 
larger that 0.1. The inset shows trace distance T> in function 
of aril'll) for different value of 9. We randomly choose 10® 
wavefunctions in a d = 10 dimensional Hilbert space. Dashed 

lines in the inset represent the curves V = ee^. 

b 


define the accuracy in terms of the trace distance T> be¬ 
tween the correct wavefunction and the weak value 
approximation tlJw,x [H], that for pure states reduces to 
V — -s/l — |('0|'0w)|2. We first give the analytical ex¬ 
pression of D in terms of the original wave function and 
then show how T> can be upper bounded by using the 
measurement outcomes. 

As shown in SI, the relation between the exact wave- 
function '0a; and the weak-value estimate ^/Jwx given in 
© can be expressed by the following relation: 


i^w,x 


'4^x 


0 - egt/j* 

N ’ 


(5) 


with eg = 2sin^(f), A/ = ^|0 - e 0 ( 0 ^)p + and 

(7^ = (|0a;P) — K0x)P- In the previous equation is 
the “variance” of the wavefunction where the average is 
defined with respect the probability density Px — | 0 xp 5 
namely (|0xP) = Ex 10*1^ ^nd (0a;) = Ex0*l0®P- 
inserting ^ into the trace distance V we obtain: 

1) = ^, (6) 

expressing T> in terms of the original wavefunction 0a: and 
the interaction parameter 9. The previous expression 
indicates when the weak-measurement method can be 
efficiently used: indeed, when 


V<^1, (7) 

the approximate wavefunction 0m,x correctly estimates 
the wavefunction 0a;. Since eq. i§ can be inverted into 
^ e0(0x)|, for small V condition 0 

is equivalent to <C 1 (see SI for the detailed 

calculation). 

Condition 0 , however, cannot be used if the exact 
wavefunction 0a; is unknown. For this reason, we now 
present a sufficient condition for the application of DWT 


method that is expressed in term of the measured prob¬ 
abilities. As shown in SI, when the follow inequality is 
satisfied 


_ pC)) > 0, (8) 


the systematic error is bounded by T* < 9/2 (for small 
9). We note that eq. 0 is equivalent to 0m > 0 when 
the global phase of 0m,x is fixed by eq. 0. 

If condition 0 is not satisfied the DWT method is 
not guaranteed to work and a lower 9 should be choosen 
to achieve condition 0. Since 0m can be expressed in 

^ T2 _ 

term of the original wavefunction as 0m = j fo^ 

wavefunction with 0 0; 0 it is possible to lower 9 suclpthat 
condition 0 is satisfied. The wavefunctions with 0 = 0 
corresponds to the set of “pathological” wavefunctions 
for which the DWT and the DST methods can never 
be applied. Indeed, if 0 = 0 the systematic error 0 
can be easily evaluated to he V — ■sj (|0xP) that is 
independent of 9: by changing the interaction parameter 
the error cannot be lowered for such wavefunctions [1^. 
Also for DST, the proportionality constant in 0 

diverges if 0 = 0. In such case, a different momentum 
state for post-selection different from |po) must be used. 

To better evaluate the accuracy of the DWT we have 
randomly chosen 10® wavefunctions in a d = 10 di¬ 
mensional Hilbert space according to the Haar mea¬ 
sure. We calculated for diherent values of 9 the prob¬ 
ability Pm Icfi violate the sufficient condition, namely 
Pm = Prob(0m < 0). We also calculated the proba¬ 
bility px> of having an error D, evaluated by 0. larger 
that 0.1. In Figure we show the probabilities pm and 
pp in function of 9. In the inset we also show the sys¬ 
tematic error V in function of (Jp/0 for different values 
of 9. Since the distribution of A/” is peaked around 0 for 
9 < 0.5, it is possible to approximate V ~ eg^: indeed, 
dashed lines in the inset of Fig. [^represent the curves 

V — eg^. The figure shows that for low values of 9, 

w 

the DWT method fails with low probability and the sys¬ 
tematic error is limited. Indeed, if we choose 9 < 0.2 for 
the d = 10 case, we have pm < 1.75% and pp < 0.57%. 
Then, as expected, low values of the interaction parame¬ 
ter 9 are suitable for the correct application of the DWT 
method. However, as we will show in the following, such 
low 9 values lead to a larger statistical error (i.e. lower 
precision) compared to the strong measurement method. 

Precision of the DWT - An important performance 
parameter is the precision of the method, namely the 
statistical errors on the estimated wavefunction. In par¬ 
ticular, it is important to evaluate the scaling of such 
errors with the number of measurements. To this pur¬ 
pose, we evaluated the mean square statistical error 5tf 
of the DWT and DST methods, obtained by summing 
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FIG. 3. Ratio of statistical errors function of 'ipw- 

Shaded area represent the points in which the DWT is con¬ 
venient with respect to the DST method, corresponding to 
'ipw > 0 and Sipw < Sips- 
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FIG. 4. Ratio of statistical errors ., in function of T). 
Shaded area represent the wavefuntions for which the sta¬ 
tistical error of the DWT is lower than the DST method. 


the squares of the statistical error on the different ip^'- 




(9) 


As shown in SI, the ratio between the statistical errors 
5ips and 5ipw^ respectively corresponding to the strong 
and weak method, can be approximately bounded by: 


6ipw 


< 


sm Uq 



(2d - 5)^2 + 2^ + 8 - 2/d 


(2d - 1)^^ 


260(1 -^- 2 ^ 2 ) 

( 10 ) 

where Oq is the interaction parameter used for the weak 


measurement. The terms sin Oq in eq. (10) shows that low 


values of Oq correspond to a lower precision (i.e. larger 
statistical errors) of the DWT with respect to the DST 
method. In the statistical analysis, we compared the two 
method by hxing the number of repetition N of the ex¬ 
periment: in the DWT or DST method, N/2 or N/3 
repetitions are used for each basis respectively. This is 
the origin of the ■\f^j2 factor in eq. (10). 

For a complete demonstration of such feature we calcu¬ 
lated the exact ratio for 10® randomly chosen wave- 
Unctions and compared it with the success parameter 
ipw and the systematic error V. The results are shown 
in Fig. and Figure show that, when the suffi- 



0 


{V) 


FIG. 5. Mean values of and T> averaged over 10® 

random wavefunctions in function of 9. 


ciency condition for applying the DWT is satished, (i.e. 

> 0), the statistical errors of the DWT are typi¬ 
cally greater then the errors of the DST. An approximate 
trent of the ratio Sips/Sipw can be obtained by noticing 
that, since J\f ^ ip, we can approximate ipw ~ "0 — ce/'^- 
Dashed curves in Fig. [^represent the r.h.s. of eq. (10), 

with ip replaced by | {ipw + repro¬ 

duce the behavior of the ratio Sips/Sipw 

To further prove that the DST precision is typically 
greater than the DWT one, we plot in Fig. [^the same 
ratio Sips/Sipw in function of the exact trace distance 
V: for low systematic error V, the statistical errors of 
the DWT are typically greater then the errors of DST. 
Equivalently, statistical errors of the DWT are reduced 
only as the systematic errors increase. Fig. shows that 
the DST precision overcomes the DWT one in most of 
the cases in which the DWT is accurate. 

To better appreciate the above results, we plot in Fig. 
5 the mean values of and T> averaged over 10® ran- 
aom wavefunctions in function of 9. The plot in Fig. 
shows again that in order to lower the trace distance T> 
it is necessary to decrease 9. However, decreasing 0, the 
statistical error Sipw becomes larger than dips- 

Mixed states - The DWT can be generalized to de¬ 
termine the density matrix p of mixed states, as shown 
in [12]. To directly measure p the same method de¬ 
scribed for pure state can be used, with the extra re¬ 
quirement that the strong measurement on momentum 
should be performed in all the momentum states \p) — 
x), while the pointer is measured is the 
|±)-p, \R)'Pi \L)v states (as done for the pure state \ip)x)- 
We indicate by p^ the density matrix that is recon¬ 
structed by the DWT and that approximates the correct 
matrix p. As shown in SI, it can be expressed as 


p'^ = 


cos 9 


[p + {cos 9 - l)D] , 


( 11 ) 


with D a diagonal matrix whose element are equal to the 
diagonal of p, namely Dx,y — dx,yPx,x- By evaluating the 
accuracy of the DWT in terms of the trace distance T) 
between p and p^ we obtained 


2cos0 




( 12 ) 
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Also in this case, the larger is the larger is T> and 
the lower is the accnracy in the estimation of p by the 
DWT. Similarly to what we have shown for pnre states, 
by performing an extra measnrement of the pointer in 
the |l)p state, the exact expression of the density matrix 
can be obtained for any valne of 9 also in the case of 
mixed states (see SI). 

Conclusions - We have demonstrated that, in order 
to achieve a direct measnrement of the wavefnnction, 
weak measnrements are not necessary. Indeed, we have 
shown that by nsing strong measnrements, in which a 
large entanglement is achieved between the system and 
the pointer, a better estimation of the wavefnnction, in 
terms of precision and accnracy, can be obtained for ran¬ 
dom matrices in most cases. Onr method allowed ns to 
derive a snfRcient condition for the applicability of the 
Direct-Weak-Tomography. We believe that onr resnlts 
give a deeper nnderstanding of the meaning of the weak- 
valne for the estimation of the wavefnnction. 

We thank P. Villoresi of the University of Padova and 
L. Maccone of the University of Pavia for nsefnl dis- 
cnssions. Onr work was snpported by the Progetto di 
Ateneo PRAT 2013 (CPDA138592) of the University of 
Padova. G.V. also acknowledge the Strategic-Research- 
Project QUINTET of the Department of Information En¬ 
gineering, University of Padova. 
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Supplementary information: 

Strong measurements give a better direct measurement of the quantum wave function 


Derivation of the wavefunction by generic strength measurement 


We here demonstrate the relation given in eq. (4) of the main text. Let’s consider a generic interaction parameter 
6 and the input state |4^in) = l'0)x ® |0)-p- In this case the unitary operators Ux{0) becomes 


Ux{0) — tx ® [(1 “ cos0)l7r + ^sin6*ay] , 

and the (unnormalized) pointer state after the momentum post-selection is given by 


(SI) 


\^)r = (Po|t/x(^)|^in) = ^[^1^) + , 


(S2) 


were we have dehned the (unnormalized) state |x) = (cos0 — l)|iL) + sin0|T/) and ^ ^'Ex tl^x- As indicated in the 
main text, it is possible to choose the phase of the wave function such that t/j — By dehning eg — 2sin^ the 
probabilities of measuring in the different pointer states are given by 


p(x) 

~dl 


1 


Pi sin 


- 2e0'ijj^e{'ijjx) + 

0 


~d 


pix) _I 

+ d 

1_ 

- {eg - sm 9 )ifj^e{'i/jx) + (1 - sin6»)e0|^x|^ 

~ d 

^ + 0 5 Re(^x) 

p{x) 

d 

1 

J 1_1 

- {eg + sin 6 »)^ 5 Re(^x) + (1 + sin6»)e0|^x|^ 

~ d 

1 1-1 

H 

q; 

1 

J 1_1 


(S3) 


pix) _£ 

~d 


pix) _£ 

~d 


+ 0 tlj Qm{tfjx) + eg ^ 5 Re(^a;) 


^ - sin^^ Am(^a,) + eg -i/j^e{tjjx) 


^ + 0^m{'i{jx) 


^ - 0^m{ifjx) 


For low the approximate results of the r.h.s. holds (at the hrst order in 9). 

We note that, by dehning aj — (ej|0) and /3j — (ej|l), the above probabilities can be obtained by the following 
POVM: 


Pj = Tr 




(S4) 


with 


Ej = ajbo)(po| -7j|Po)(a^| 


(S5) 


and 7 j = [(1 “ cos 9)aj — sinO/dj]. 

From the previous equations (S3), using the exact results, it is possible to prove that: 


^e{t/jx) 

Am(^a,) 


d 


(x) 


2'0 sin 9 


d 


2-0 sin 9 


[P 

[p(x) 


p!*^+2tan(^)Pi^“^)], 
P^R^] • 


(S6) 


Strong measurements correspond to 6* = 7r/2. By measuring the pointer in the |ei) = |1), |e_|_) = |+), |e_) = |—), 
Icl) = \L) and |ei?) = \R) basis, the wave function can be thus derived. It is very important to stress that the result 
is exact, without any approximation. 
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(x) 

If we consider the weak-value approximation, then the approximate valnes of Pj can be nsed. In this case 


MM = - -P-"’) 

20tp 

20tp 


(S7) 


Relation between weak and strong value 


Let’s now derive the the relation between the correct wave fnnction and the weak-valne estimate To 

evalnate the wave fnnction it is necessary to estimate the parameters + 2tan(|)Py^^ and — 

pA) _ pA) ^j^at the wavefnnction is obtained by normalization: 




Ax + 

M 


with M. — + ^x) — ^ ■ On the other hand, the weak valne wave fnnction is given by 


I Aw,x 

H^W^x — 


iB 


W,x 


M 


w 


(S8) 


(S9) 


with the parameters given by A]4r,x — — P}_^^ — A^[3?e('0a;) “ ^w,x — — Pr’ — M.Am.{xl)x) and 

Mw — X + ^Wx)- comparing the two resnlts we can write 


)(») _ 


_ p(®) 


)(») _ 




(SIO) 


with A/* determined by the normalization of 'i[^w,x- 




|2 _ 


^2 - 2e0'ilj^e{'ilJx) + e^(l^xp) 


(Sll) 


|^- e ,(^.)|2 + e 2 a 2 . 


The average is dehned with respect the probability density dehned by the wave fnnction, px — |'0xP, namely (It^xP) = 
'^x {^e{xjjx)) — J^x ^g('0x)|'0xP- We have thns demonstrated eq. (5) of the main text. 

We now show that the trace distance between '0m,x and xjjx can be bonnded by knowing tl^w,x- The trace distance 
can be exactly evalnated if we know the correct wave fnnction tfjx, by 


V = 




(S12) 


The relation between P and can be inverted by sqnaring the previons eqnation, namely — 
By resolving for ega-ip we obtain 




=p^\x{j—eg{x{jx)\ that for low V can be approximated by ega^ ~ P|'0—e6»('0, 


|V'-ee(V'cc)L + (eeo-^)2 

x)|, 


snch that V 




. Then, for small P, condition P <C 1 is eqnivalent to 


^0 ^xp 


<C 1. 


Ib-e0(bx)| _ |■i/’-ee(p^)| 

The parameter P can be bonnded by knowing xpw = Y^x '^w,x- We note that the global phase of xl)w,x is hxed by 
(SIO). By nsing (SIO), we have xpw — — eg)/Pf. When xpw > 0 we can conclnde that xp > -s/ee allowing to bonnd 

the parameter P. Indeed, since \{xl)x)\ < 1 and eg < 1, the condition xp > ^Jeg implies \xp\ > eg\{xpx)\ and 

- ^9\{ipx)\f > {^-£9^ > e0(l - ■ 

Finally, since dp < l/\/2 we can conclnde that 

1 - 1 - \/p9 


(S13) 


xpw > 0 


P = 


< 


|■0-e0(■0x)|^ 

4-1 


+ 1 


< 


0 (1 yW)^ _|_ ^ 


e/2 


(S14) 


























where the last approximate result holds for low 0. The previous relation proves that the condition > 0 gives an 
upper bound on the systematic error V. By eq. (SIO) the sign of il)w is equal to the sign of 


then 


equation (S14) proves eq. (8) of the main text. 


Analysis of the precision of the DWT and DST methods 


It is useful to introduce the following average error Si/j, obtained by summing the squares of the statistical error on 
the different 


^4’=JpAA 


with |(5'0xP = 55Re('0x)^ + 

In general, for a wavefunction written as ijjx — with M. — \/Ylxi^x + ^x) have 


(S15) 


Sijj — 





)5Al + (1 


-^x 


)5B^x] 


Let’s now evaluate the above expression for the two methods, the DWT and the DST. 


(S16) 


DWT 


Let’s now evaluate such average error 5ijj in the weak measurement case. Let’s consider to repeat the experiment 
N times. For the weak value we need to measure in the {|I/), \R)} and the {|+), |—)} basis. Let’s suppose that N/2 
measurements are used for the hrst basis and N/2 per the remaining basis. We indicate with tilde the estimated 
parameters obtained ofter N measurements. The estimate for the polarization probabilities Pj are 


Pj- 


N/2 


N- 




(S17) 


since nj —> PjN/2 in the large N limit. From now on we indicate with a right arrow the asymptotic behavior in the 
large N limit. The variance of nj is equal to nj due to Poissonian statistic. Then 




N 


‘^Pj 

N 


(518) 

(519) 
with the factor 

M-w determined by normalization M.w — ^jYhxi^wx + Pwx)- large N limit, the estimated M.w approaches 

w —^ statistical error on the estimated Ax and Bx are given by 


^ N 

The probabilities are used to estimate the terms 

Aw,x — P+ — P- 1 Bw,x — Pl — Pr , 

from which the wave function is obtained in the large N limit as 5Re('0x) = and Am(i/jx) — ^ 

w 


to M- ^ 


SAw,x = V S^P+ + S‘^P- —^ ~ 


(S20) 


sBw,x^ys^PLPS^PF 


5Aw^3. 


Since in the large N limit we have 


given by 




5fte(^m,cc), 


AT, 


Am.{'il)w,x)^ the mean square statistical error di/w is 


5^/w = 


M 


w 


^[(1 - + (1 - 


(S21) 
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From (S20) we have 


X 


dN 


dN 


2 d' 0 ^ + 466/(1 — ip^) 


^ - 2e0^5fte(^^) + 2e0|^a;|^) 


(S22) 


By using the previous equation in (S21) we obtain, for the weak value case, 


- 1) ^2 ^ 460(1 - ^2) + 260 ^ |^M/,xp(^5Re(^x) - l^xP) 


(S23) 


where we used - \ipw,x?\'^x? > -1 and ^Y.x I'0w,xp3^e(^x) > -^. 


DST 

Let’s now evaluate the average error Sip in the strong measurement case. Again we consider to repeat the experiment 
N times. We note that in this case we need to measure the ancillary qubit in three bases, namely {\L), |i?)}, {|+), |—)} 
and {|0), |1)} basis. Then, N/3 measurements are used for each basis, such that 


6P, 


N 


(S24) 


The Ax and Bx factor in eq. (S27) are given by: 


Ax — P+ — P- + 2tan(-)Py , — Pl ~ Pr i 


(S25) 


with errors 


5A^ = 5‘^P+ + 5‘^P- + 4 tam(-) 52 p^ —^ ^2 _ 2eetp^e{tpx) + 2|^^|260(1 + 26^) 


5B, = ^5^Pl + 5^Pr 1/ ^ V - 2e0^5Re(^,) + 2e0\P^x? • 


(S26) 


In the large N limit we have ^ ^ ‘^e{'ipx)^ ^ Am{ipx), and the mean square statistical error dtp become 


<5^ = - 3f^e(^^)2)(5A2 + (1 - Am{tpx)^)6B‘^] 


(S27) 


By using (S26) we have 


J2(SAl + SBl) 

X 


Nd 


Nd 


2d'ip‘^ + 46^(1 + 60 — ' 0 ^) 


tp^ - 2e0tp^e{tpx) + 260 ( 10 ^ 1 ^) + 4:el{^e{tpxf) 


(S28) 


In the large N limit we have A4 —> ^ such that 


dtp — 


3d 


ips'mO V 4A” 


-T^\/{2d- 1)02 + 4 e 0 (i + 60 - 02 ) + 26005Re(0x) - 260 ( 10 : 02 ) - 4:e‘^{^e{ipx)^) 


(S29) 
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For the strong measurement we have d — tx j2 and = 1 such that 

&i>s = - 5)*?" + 8 - -iW) - 4(Sfe(^*0"> + 2^^ (Se(^O) 

^ ,_ (S30) 

1 //- 

< {2d - 5)t^ + 2t + & - 2/d 

since (5Re('0x)) < 1, (iV^xP) > ^ and '0 < \/d. 

The ratio between the statistical errors can be bounded by 


St/Js . . [SAf I (2d-5)02+ 20+ 8-2/d 
Sifjw y ^ ifj \] (2d — 1) 02 + 2e0(l — 0 — 20^) 


(S31) 


that is the main result (eq. (9) of the main text) due to the fact that ^ is well peaked around 1. For large 0 we have 


J\f ^ 0iy and the bound is simplihed to 


sin 


^\/I 


2<i-5 

2<i-l 


Mixed states with intermediate measure 

Let’s consider the system initially prepared in the state 

n 

p = px 0 |0)p(0|, with px ^ Px,y\x){y\ (S32) 

x,y=l 

and |0)p the initial “pointer” state. We would like to estimate the density matrix px- After the interaction with 
Ux{0) — , the state becomes 


p', = uMpUl{e) 


(S33) 


The system is then measured into the momentum state \p) — such that the remaining “pointer 




becomes: 


with elements pQo{x,p) 


Px,p = {P\p'x\p) = 


Poo{x,p) poi{x,p) 
Pw{x,p) pii{x,p) 


(S34) 


'^x,y Px,y^ 


2-Ki{y — x)p 


- 2sin2 I (px,ye 


2-Ki{y — x)p 


+ c.c.) + 4sin^ |p; 


2 Hx^x 


, Pio{x,p) 


4 sin I 


2iTi{y — x)p 


- 2sin^ |p:r,x|, Poi(ic,p) = pio(ic,p)* and pii(x,p) J sin^ dp:r,x- 


Now it possible to determine Px,y in function of the “pointer” density matrix as follows: 


U \ ^ 2,7ri(^x — y)p 

px,y oc dtan -Sx,yPn{x,p) + 2_^e ^ Pio(a;,p) 

p 

I strong measure {9 — 7r/2) 

E 27rz(x — y}p 

e d pio(x,p). 

p 


The weak value estimate is obtained at the lowest order in 9\ 

2TTi{x-y)p 


-W 

Px,y OC 


E ^7rz{x — y)p 

e ^ Pio(x,p) 


p 


(S35) 


(S36) 


(S37) 


By the above equation it is possible to express the estimated density matrix p^ in terms of the correct density p as 

1 


w 

P = 


cosd 


[p + (cosd — l)i4] 


(S38) 
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with D a diagonal matrix whose element are eqnal to the diagonal of p, namely Dx,y — Sx,yPx,x- 

To determine the terms pio in ( |S37| ) it is necessary to measnre the pointer into fonr states |+)p, |—)'p, |T)-p and 
\R)'P. Indeed, by dehning = {i\p''\i)^ we have the DWT relations: 


Piq{x,p) 


1 - 
2 - 




-P 


P,p) 


(S39) 


The DST allows to determine the exact density matrix p by also measnring the pointer into the state |l)p. Indeed, 
to determine p by eq. (S36) it is necessary to evalnate also the the terms pi\. It is easy to show that snch terms can 


be evalnated by measnring the pointer into the states ll)^: 


pn(x,p) = Pp> = P’ 


(S40) 


To snmmarize, for mixed states the procednre is similar to the one performed with pnre state. The initial state 


(S32) is evolved according to the interaction Ux{0) between the system and the pointer state. The system state is 
strongly measnred in a given momentnm state \p) snch that the pointer is left into a mixed two-level system p^p given 
by eq. (|S35 ). By performing a standard tomography on the pointer, namely by projecting it into 1+)^, |—)p, \L)p 


and |i?)p the pointer state can be obtained. Then, the exact initial density matrix can be derived by eq. (1 ^ . 













